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Abstract. It is clear from stellar evolution and from observations of WR stars that massive stars are releasing metal-enriched 
gas through their stellar winds in the Wolf-Rayet phase. Although Hii region spectra serve as diagnostics to determine the 
present-day chemical composition of the interstellar medium, it is far from being understood to what extent the Hii gas is 
already contaminated by chemically processed stellar wind. There fore, we analyzed our models of radiative and wind bubbles 
of an isolated 85Mo star with solar metallicity iKroser et all2006h with respect to the chemical enrichment of the circumstellar 
Hii region. Plausibly, the hot stellar wind bubble (SWB) is enriched with ^^N during the WN phase and even much higher with 
^^C and ^^O during the WC phase of the star. During the short period that the 85Mo star spends in the WC stage enriched SWB 
material mixes with warm Hii gas of solar abundances and thus enhances the metallicity in the Hii region. However, at the end 
of the stellar lifetime the mass ratios of the traced elements N and O in the warm ionized gas are insignificantly higher than 
solar, whereas an enrichment of 22% above solar is found for C. Important issues from the presented study comprise a steeper 
radial gradient of C than O and a decreasing eff'ect of self-enrichment for metal-poor galaxies. 
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1. Introduction 

Hii regions are considered as the most reliable targets for the 
determination of the present-day abundances in the Interstellar 
Medium (ISM). Nevertheless, heavy elements are released by 
massive stars into their surrounding Hii regions. The knowl- 
edge of its amount is of particula r interest for the obser- 
vation of very metal-poor galaxies. iKunth & SargentI {1986) 
e.g. discussed the problem of determining the heavy-element 
abundance of very metal-poor blue compact dwarf galaxies 
from emission lines of Hii regions in the light of local self- 
enrichment by massive stars. This can happen during two 
stages in their evolutionary course: In the Wolf-Rayet (WR) 
stage, when the stellar wind has peeled oflT the outermost stel- 
lar layers, and by supernovae of type II (SNell). Because of 
their energetics, stellar winds and SNell contribute their gas to 
the hot phase. Therefore, it is of high relevance for observa- 
tions and the addressed question of the presented exploration 
to study the amount to what it can be mixed into the diagnosed 
Hii regions. 

As discussed by IChiosi & MaedeJ (Il986l) . massive stars 
with an initial mass greater than 25 M© evolve into WR stars 
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with enhanced mass-loss rates and chemically enriched stellar 
winds. 

lAbbott & Contil (Il987l) concluded from stellar evolution 
models that WR winds strongly contribute to the Galactic en- 
richment of ^He, ^^C, ^^O, and ^^Ne, while contributing only 
moderately to the enrichment of ^"^N, ^^Mg, ^^Mg, and ^^O. 

For metallicities less than solar two effects reduce the heavy 
element release by WR stars: First, the lower the metallicity the 
more massive a star has to be to evolve through the WR stages. 
Therefore, t he number of WR sta rs decreases with decreasing 
metallicity. ISchaller et all (ll992L hereafter SSMM) found that 
with metallicity Z = 0.001 the minimal initial H-ZAMS mass 
for a WR star is > 80 M©. At Z = 0.02 the minimal initial 
mass is > 25 M©. Second, the lower the metallicity the shorter 
are the WR lifetimes, and not all WR stages are reached. At 
solar metallicity WR stars enter all three WR stages (WNL, 
WNE, WC), whereas at Z = 0.001 only the WNL phase is 
reached ISSMM) . The WR lifetime of an 85 M© star, e.g., is 
%R = 0.2 04 X 10^ yr at Z = 0.001 and %r = 4.008 x 10^ yr at 
Z = 0.02 (ISSMMh . 

That WR stars play an important role for car bon enrich- 
ment o f the ISM at solar metallicity was shown bv lPrav et al.l 
(I2OO3I) . Their models predict that the C enrichment by WR stars 
is at least comparable to that by AGB stars while the enrich- 
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ment of N is dominated by AGB stars and the enrichm ent is 
dominated by SNell. At other metallicities lPrav & Tou3 (Eo03) 
found an increase in carbon enrichment with increasing metal- 
licity but a decrease in oxygen enrichment. The N enrichment 
by WR stars is negligible compared to that by other sources. 
Additionally, the total amount of N ejected by a WR star gen- 
erally decreases with decreasing metallicity. 

iHirschi etal .1(12004^)051 showed that the consideration of 
stellar rotation in the Geneva stellar evolution code increases 
the yields for heavy elements by a factor of 1.5 - 2.5 for 
stars between 15 and 30 M©. For the more massive stars ro- 
tation raises the yields of ^He and other H-burning products 
like ^"^N, but the yields of He-burning products like ^^C are 
smaller. Additionally, for stars with M > 60 M©, the evolution 
differs from that of non-rotating stars by the following manner: 
Rotating stars enter the WR regime already in the course of 
their main- sequence. 

As an extent to stellar evolution models mentioned above 
hydrodynamical simulations of radiation and wind-driven Hii 
regions can trace the spatial distribution of the elements ejected 
by the star, and, therefore, we can study the mixing of the hot 
chemically enriched stellar wind with the warm photoionized 
gas. 



2. The model 

We analyze the model of an 85Mq star in a series of radiation 
and wind-dr i ven Hii regions arou nd massive stars (iFrever et alJ 
I2003L 12004 lKrogeret alJl2006[ her eafter: Paper I, Paper II, 
Paper III) as described in Paper IIlL A descripti on of the nu- 
merical method and further references are given in lPaper J . The 
time-dependent parameters of the 85 Mq star with "standard" 
mass-loss and solar metallicity (Z = 0.0 2) during its H-MS and 
its subsequent evolution ar e taken f rom SSMM, who used the 
relative abundance ratios of lAnders & Grevesse, (.1989.) for the 
heavy elements. 

An 85 M© star with solar abundance enters its WN stage 
at an age of ^ = 2.83 Myr. The WR star enriches the 
combined SWB/Hii region with ^^C, ^"^N, and ^^O. In our 
2D simulations we study the spatial distribution of these 
tracer elements inside the SWB/Hii region as well as their 
time-dependent mass ratios as Mi2(^/Mtotai, Mi4|y[/Mtotab and 
Mi6Q/Mtotai in two different temperature regimes: "warm" gas 
(6.0 X 10^ K < r < 5.0 X 10^ k) accounts for the gas of the Hii 

region, and "hot" gas, whose temperatures (r > 5 x 10"^ are 
reached inside the SWB. Additionally, for the Hii region diag- 
nostic a degree of ionization of > 0.95 must be reached. The 
undisturbed ambient medium outside the Hii region as well 
as the material swept-up by the expanding Hii region are not 
counted since they are not ionized. The observable abundances 
of Hii regions are derived for the warm phase only. 

The mass fractions of ^^C, ^^N, and ^^O in the circumstel- 



lar gas as well as o f the star itself are initially set to the solar 
values, taken from lAnders & Grevessd (Il989h . 4.466 x 10"^, 
1.397 X 10"^ and 1.061 x 10"^, respectively, in relation to H. 



3. Results 

The spatial distribution of the tracer elements is shown as 
log(p (element) /p (element, solar)) at different stages of stel- 
lar evolution, where "element" means ^^C (Fig.[i] middle col- 
umn), ^^N (Fig.[ll left column), or ^^O (Fig.[l] right column), 
respectively. All plots cover the whole computational domain 
of 60 pc X 60 pc. The gas density is overlaid as a contour plot. 
Before the onset of the WN phase, these density ratios are unity 
by definition. With the onset of the WN phase the combined 
SWB/Hii region is chemically enriched by the WR wind. 

We study the temporal evolution of ^^C, ^"^N, and ^^O, re- 
spectively, normalized to solar values in the two different tem- 
perature regimes of the combined SWB/Hii region: For the 
hot gas of the SWB the mass ratios are depicted in Figure |2l 
Figure [3l shows it for the "warm" Hii gas (see §|2|). 

3.1. Chemical enrichment of the stellar wind bubble 

In the hot gas of the SWB the abundance of ^"^N at t = 2.99 Myr 
(close to the end of the WN phase, Fig.[T] upper left plot) that 
lasts only for ^ 0.17 Myr, reaches values in the range of 0.1 < 
log(^4N/solar) < 0.98. is distributed in the whole SWB 
except in the Hii material enveloped by the SWB. ^^C and ^^O 
can be found at the same places, but their concentration is as 
low as log(^^C/solar) ^ -1.4 and log (^^0/ solar) ^ -1.5. 

During its WN phase the star releases 0.143 M© of ^"^N, 
which is more than half of its total ^"^N release, but almost no 
extra ^^C or ^^O is supplied. 

With the onset of the WC phase at t = 3.005 Myr the 
^^N release ends and the main enrichment of the SWB with 
^^C and ^^O starts. Due to the beginning WC phase with its 
high mass-loss rates and high terminal wind velocity the pre- 
viously ejected material in the SWB including the emitted 
^"^N is compressed into a shell-like structure. Therefore, after 
t = 3.005 Myr ^^C and ^^O are distributed in the whole SWB 
and ^"^N can mostly be found in an expanding shell (see lower 
left plot of Fig.[l] at ^ = 3.22 Myr). The density contour line at 
a radius of 28 pc in the upper middle and upper right plots of 
Fig.[T](^ = 3.06 Myr) marks the position of the reverse shock 
which borders the free-flowing wind zone. At 3.06 Myr the ra- 
dius of the reverse shock has reached its biggest extent. 

During this evolutionary phase (see upper middle and upper 
right plots in Fig.[T] at 3.06 Myr) the star ejects large amounts 
of ^^C and ^^O: The log (^^C/ solar) value reaches ^ 1.86 in- 
side the highly enriched SWB and the log (^^0/ solar) maxi- 
mum value is ^ 1.35. 

At the end of its lifetime at ^ = 3.22 Myr the 85 M© star 
has supplied 0.28 M© of ^^N, 13.76 M© of ^^C, and 11.12 M© 
of ^^O, which are contained in the combined SWB/Hii region. 
This situation is shown in the lower plots of Fig. [I] The reverse 
shock is now located at a radius of 9 pc, the hot gas of the SWB 
fills the area up to a radius of 39 pc. Almost all of the ejected 
i^C, i^N, and ^^O is located in the SWB's volume. 

In the outer zones of the SWB the highest abundance of 
^^C (see lower middle plot of Fig. [Q is locally reached with 
^ 1.85. For radii r < 15 pc the value is ^ 1.74. 
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model: 605 time: 2.99E + 06 yr 
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Fig. 1. Concentration of ^"^N (left plots), ^^C (middle), and ^^O (right plots). Upper row: t = 2.99 Myr (top left), near the end of 
the WN phase, and at ^ = 3.06 Myr (middle and right), after the onset of the WC phase. Lower row: At the end of the stellar 
lifetime (t = 3.22 Myr). 



The log (^"^N/solar^ ratio (Fig.[T] lower left) reaches ^ 0.55 
in the outer zones of the SWB, where the material ejected dur- 
ing the WN phase is located. In the inner parts of the SWB only 
^^C and can be found. 

The maximum of the log ^^^0/solar^ value (Fig.[n lower 
right) of ^ 1.65 is reached at radii r < 15 pc. Within the outer 
parts of the SWB the values range from 0.5 to 1.4, highly de- 
pending on the location, since the abundance within the em- 
bedded clumps differs from that of the hot SWB material. 

The temporal evolution of the ^^C abundance in the hot gas 
phase is plotted in Fig.|2las dash-dotted line. The plot starts at 
t = 2.83 Myr, when the WN phases are reached. Additionally, 
the onset of the WC phase at ^ = 3.005 Myr is i ndicated. During 
the WN phase the stellar parameters (IsSMMI) specify the ^^C 
abundance of the WR wind to about 1/14 solar. Due to the mix- 
ing of the WR wind with the main- sequence wind the averaged 
mass ratio in the hot gas of the SWB is only 0.7 times solar 
during the WN and early WC phase. The ^^C enrichment starts 
with the onset of the WC phase. After the newly ejected mate- 
rial passed the reverse shock the ^^C abundance in the hot gas 
reaches 2.1 times the solar value at t = 3.02 Myr. This value 
increases further to 38 times solar at the end of the stellar life- 
time. 

The ^^N abundance is depicted in Fig. |2| as solid line. It 
rises as the enriched free-floating stellar wind passes the re- 
verse shock, where the material is heated to several 10^ K. 
Since in the WC stage the stellar wind contains no ^^N, the ^^N 



abundance in the hot SWB decreases strongly after it reached 
its maximum value of 3.4 times solar at ^ = 3.0 Myr. At the end 
of the stellar lifetime the ^"^N abundance in the hot gas phase is 
only 0.8 times solar. 

Because the enrichment with ^^O starts together with ^^C 
at the onset of the WC phase and due to the high mass-loss rate 
in the WC phase the 0-abundance in the hot gas phase rises 
steeply after the material passed the reverse shock. After t = 
3.02 Myr the ^^O mass ratio in the hot gas becomes supersolar 
and increases until the end of the stellar lifetime up to 15.7 
times solar. 

3.2. Chemical enrichment of the Hii region 

During the WR phases of the 85Mq star the SWB expands into 
the surrounding Hii region without being bordered by a shock 
front. Therefore, the hot SWB gas structures the Hii region and 
encloses Hii gas which is thus compressed into clumps. This 
embedded Hii material contributes to the warm gas phase as 
well as the Hii layer around the SWB. 

At the end of the stellar lifetime the ^^C abundance (lower 
middle plot in Fig.[i]) in the Hii layer reaches a local maximum 
of log(i2c/solar) ^ 0.17. It contains a small fraction of ^"^N, 
the log (^"^N/ solar) values there are up to 0.01 (Fig. [T] lower 

left). The log (^^0/solar) values (Fig.[l] lower right) range from 
Oto 0.03. 
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Fig. 2. Time-dependent abundance of ^^C, ^"^N, and ^^O in the 
hot gas phase. For details see text. 
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Fig. 3. Same as Fig.|2l for the warm Hii gas phase. For details 
see text. 



The temporal evolution of the ^^C abundance is plotted as 
dash-dotted line in Fig.|3l As in Fig.|2l the plot starts with the 
beginning of the WN phases, the vertical dotted line indicates 
the onset of the WC phase. In the warm Hii gas the averaged 
^^C abundance reaches supersolar values after t = 3.02 Myr, 
thus the enriched material needs about 1.5 x 10"^ yr for its way 
from the stellar surface, through the SWB, until it reaches the 
Hii region and mixes with the warm Hii gas. It rises steeply 
until reaching a value of 1.21 times solar at ^ = 3.15 Myr and 
is 1.22 times solar at ^ = 3.22 Myr. 

The ^^N mass fraction (solid curve in Fig. O in the warm 
gas phase rises to supersolar values after ^ = 2.88 Myr when the 
^^N has passed the SWB and cools down to Hii temperatures. 
As discussed for Fig.[ll after ^ = 3.0 Myr the ^"^N-enriched ma- 
terial can only be found in the shell-like structure in the outer 
zones of the SWB. The ^^N abundance in the warm Hii gas 
reaches its maximum of 1.008 times solar at ^ = 3.11 Myr and 
has a value of 1 .006 times solar at the end of the stellar lifetime. 

For the warm gas the ^^O mass ratio is given as dashed 
line in Fig. [3l Like ^^C it reaches supersolar values after 



t = 3.02 Myr. Until the end of the stellar lifetime it in- 
creases slightly but reaches only 1.06 times the solar value at 
r = 3.22 Myr. 

4. Discussion 

This chemical analysis of the wind and radiation-driven Hii re- 
gion around an 85 M© star with solar metallicity provides first 
quantitative conclusions to what extent the C, N, 0-enriched 
WR winds contribute to the observable abundances of the sur- 
rounding Hii region. 

When comparing the C, N, O abundances in our simula- 
tion with observed values, one has to take into account that the 
overall lifetime as well as the WR lifetime of the 85 Mq star 
are extremely short compared to le ss massi ve stars. Also, the 
choice of the stellar parameters of ISSMmI influences our re- 
sults, since other evolutionary tracks would provi de other W R 
lifetimes as well as other mass-loss rates (see, e.g.. lPaper IIll) . 

At the end of the stellar lifetime, the ^^C abundance in the 
warm gas phase amounts to 22.3% above solar, whereas we 
find only 0.6% supersolar for ^^N, and 5.5% supersolar for ^^O. 
These values measured in the warm gas phase are the quantities 
which should be compared with observed Hii regions' emission 
spectra. 

On the other hand, the hot gas of the WR wind bubble is 
highly enriched with ^^C and ^^O since the onset of the WC 
phase, while it was significantly enriched with ^"^N during the 
preceeding WN phase. 

From these models we conclude that the enrichment of the 
circumstellar environment with ^"^N and ^^O by WR stars might 
be negligible if the 85 Mq star is representative for massive 
stars passing the WR stage. Only for ^^C the enrichment of the 
Hii gas is significant. 

Since the occurrence of a WR phase is strongly metal de- 
pendent the enrichment with C should also depend on metal- 
licity like e.g. on O. One should therefore expect that the C 
gradient of Hii regiones in galactic discs is steeper than that 
of O. And indeed, Esteban et al. (2005) found A log (C/0) = 
-0.058 ±0.018 dex kpc"^ for the Galactic disk. 

In metal-poor galaxies one would expect less chemical self- 
enrichment because the stellar mass range of the WR occur- 
rence is shrinked and shifted towards higher masses and the 
WR phases are shorter. 
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